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Abstract: It has been recently found that the heavy quark-antiquark QQ pair multiplicity, 
in certain phase space region (QQ at short distance, soft and with small velocity), satisfies 
an evolution equation formally similar to the BFKL equation for the high energy scattering 
amplitude. We find the exact solution of the QQ-equation and discuss the differences with 
the BFKL scattering amplitude 
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1. Introduction 

A way to reveal properties of QCD radiation consists in resumming logarithmically en- 
hanced terms in the perturbative expansion for observables at short distances [1] (loga- 
rithms in Feynman diagrams come when vertices are in collinear or infrared configurations) . 

Recently attention has been attracted [2] on short distance observables for jet physics 
which have only infrared logarithms, to leading order. An observable of this type, consid- 
ered in [2], is the distribution in E ou t, the energy deposited by QCD radiation into a region 
away from all jets. Here one resums (to leading order) powers in a s lnQ/E out with a s the 
QCD coupling and Q the hard momentum scale (e.g. the center of mass energy in e + e~ 
annihilation). Another observable, considered in [3], is the multiplicity of a heavy quark- 
antiquark system QQ of mass M in a phase space region in which collinear singularities 
are absent. Here one resums (to leading order) powers in a s hiQ/ M.. These logarithms 
originate, in the parton language, from successive gluon branching 1 and are strictly related 
to the non-Abelian structure of QCD vertices. The resummation of these "non-Abelian" 
logarithms raises various surprises. 

The first surprise is that standard collinear factorization methods [1] do not work 
here [2] and these non-Abelian logarithms cannot be resummed simply into Sudakov fac- 
tors. The way used to resum them was to use the distributions for multi-soft gluons emitted 

1 These type of logarithmic terms are present also in all non-global jet-shape observables which have 
also collinear singularities such as the Sterman- Weinberg distribution [4]. Here their contribution is beyond 
leading logarithmic order. 
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in any angular configuration. They are known [5], but only in the planar approximation 
in SU(N C ). In the case of the -E ou t-distribution, the resummation is performed by numer- 
ical methods [2] or by a non-linear evolution equation [6] which reflects the non-Abelian 
structure of the QCD branching. In the case of the QQ-multiplicity, the resummation is 
performed [3] by a linear equation. 

The second surprise is that the equations resumming these non-Abelian logarithms 
are similar to the equations obtained in high energy physics in which one resums leading 
logarithms in the energy (rapidity). In particular, the equation for the QQ-multiplicity 
is similar to the Balitski-Fadin-Kuraev-Lipatov (BFKL) equation [7] for the high-energy 
elastic scattering amplitude. However, this similarity is only formal since the relevant 
momentum configurations of the soft gluon ensemble are completely different. For the QQ- 
multiplicity all gluon angles, 9ij, are of comparable order while their transverse momenta, 
qu, are strongly ordered. In the high energy case the opposite is true: comparable qu and 
strongly ordered 6ij. In [3] one finds the discussion of this point. 

The differences in the relevant phase space for of the soft gluon ensemble have im- 
portant consequences. First of all the QCD coupling enters differently the two equations. 
This is due to the fact that the arguments of the QCD running coupling for the soft gluons 
are given by qu which are ordered (QQ case) or comparable (BFKL case). As a result 
the expansion parameter for the leading order distributions is r ~ a s (Q)\nQ/AA in the 
QQ case and r ~ a s Y in the BFKL case, with a s the fixed coupling and Y the rapidity 
(logarithm of energy). This implies that the asymptotic regime for large r is physically 
accessible for the high energy case while it is quite extreme for the QQ case. 

A second difference is that the variables of comparable order are dimensionful {qu in 
the BFKL case) or dimensionless (9ij for the QQ case). Therefore, while qu are unbounded 
(total energy much larger than qu), angles are obviously bounded. As a result the BFKL 
equation is equivalent to a diffusion equation in a translational invariant system (using 
the logarithm of impact parameter), while for the QQ-ecjviation there is a boundary and 
translational invariance is lost. Only in infinitely boosted frames one could try to neglect 
the angular limitation (see [3] ) so that the two equations would become exactly the same. In 
this case one could use all results of high energy studies [8] to describe the Q(5 _mu ltiplicity. 

In this paper we directly study and explicitly solve the QQ-multiplicity equation. We 
do not need to use a boosted frame, we work in a general frame including the e + e~ center 
of mass. Moreover the exact solution is valid for any value of the expansion parameter 
r ~ a s (Q) In Q / M not only in the asymptotic regime. We compare this solution with the 
BFKL high energy amplitude. We find that the effect of the finite range in 6ij is relevant 
even at large r. Here both distributions are dominated by the Pomeron intercept but 
have different prefactors. It could be considered as a third surprise that the mathematical 
structure of this equation is quite rich and deserves a special study [9] . 

In section 2 we describe the observable (QQ-multiplicity) and report the evolution 
equation in [3] resumming leading logarithms. In section 3 we discuss the difference between 
the QQ-multiplicity and BFKL equations. With the aim to learn the strategy to discuss 
the solution of the QQ-equation, in section 4 we introduce and study a simple solvable 
model which bears features which are similar to the QQ-multiplicity equation. In section 
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5 we present the exact solution of the QQ-equation (details are given in Appendix A). In 
section 6 we summarize the properties of the solution. 

2. The observable and the equation 

The observable introduced in [3] is the multiplicity distribution of heavy quark-antiquark 
system QQ of mass A4 and momentum \k\ produced in e + e _ annihilation with center of 
mass energy Q (similar observable could equally well be defined for deep inelastic scattering 
or hard events in hadron-hadron collisions). The distribution was studied for Q^>M where 
the perturbative coefficients are enhanced by powers of InQ/M. Moreover the QQ system 
was considered near threshold (A4 2 — 4M 2 ~M 2 with M the heavy quark mass) and with 
small velocity v = \k\/E k so that there are no collinear singularities. As a consequence, 
the leading logarithmic contributions (a 2+n ln n Q/M) are obtained by considering soft 
secondary gluons emitted in all angular region off the primary quark-antiquark. The QQ 
system results from the decay of one of these soft gluons (actually the softest one). 

As shown in [3], to leading logarithmic order, the QQ-multiplicity distribution factor- 
izes into the inclusive distribution / for the emission of the soft off-shell gluon of mass A4 
and momentum \k\ and the distribution for its successive decay into the QQ system 



E k dN a 2 s C F jM 2 - 4M 2 M 2 + 2M 2 _ J 



dM 2 d\k\ 3vr 2 M 2 V M 2 M 2 

with / a function of a s InQ/M. and of v (we work in the e + e~ center of mass). Instead of 
a s lnQ/M, it is convenient to introduce the equivalent variable 2 

[Q dq t N c a s (q t ) _ 2N 2 f InQ/A \ 

' J M Qt * " nN c -2n f \lnM/Aj ' 1 " ) 

with N c the number of colours, nj the number of quark flavours and A the QCD scale. 

We consider the e + e~ annihilation in the center of mass so that, for soft secondary 
radiation, the primary quark-antiquark are back-to-back. However, to resum the leading 
logarithmic contributions (via use of recurrence relations) , one needs to consider a general 
case in which soft gluons are emitted off a hard dipole forming an angle 9. The inclusive 
distribution for a soft off-shell gluon of mass A4 and velocity v emitted off a dipole forming 
an angle 6 will be denoted by 

I(p,T,v), p=i(l-cosfl). (2.3) 

The physical case corresponds to = ir. For Q~A4 the inclusive distribution is the Born 
contribution given by (r = 0) 

il n ^ f dQk v 2 (l-cos9 ab ) 

l{p,0,v)= / -— — — — — = pf{p,v , 2.4) 

J 47r (1 - vcos6 ak )(l - vcos9 kb ) 

with p = 5(1 — cos6* a b). For M<Q secondary radiation starts to contribute. 
2 To reconstruct the argument of a s one needs to go beyond leading order in the soft limit. 
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For v = l there is a collinear logarithmic divergence in (2.4) which needs to be regular- 
ized and resummation of secondary radiation leads to the standard multiplicity function [1], 
a distribution with both collinear and soft logarithms. Here we have v < 1 and only soft 
logarithms are involved. 

To resum the secondary radiation to leading order one needs the multi-soft gluon 
emission distribution in all angular configurations as given in [5]. This distribution can 
be represented as a successive branching in which a general dipole of massless momenta qi 
and qj emits a soft massless gluon q 



dwij(q) 



N c a s (q t )dq t dfl q 2 p^ 



dr • 



Pij 



(2.5) 



nq t 47T PigPgj 2lT PigPgj 

From this branching distribution, taking into account (strong) energy ordering, one de- 
duces (see [3]) the following evolution equation resumming soft logarithms (neglect r and 
v dependence) 

d&>a Pab 



d T I(p ab ,T) 



I 



■[I(Paq,r) +I(p qb ,T) - I(p ab , t) 



(2.6) 



2VT PaqPqb 

with pij = ^(1 — cos 6{j). The v dependence, here understood, is coming from the initial 
condition (2.4) at r = 0. The first two terms in the integral originate from emission out of 
aq and qb dipoles which result from the branching of the parent dipole ab. The third term 
originates from virtual corrections in Feynman diagrams. 

Since the off-shell gluon decaying into the QQ system is the softest one, the initial 
condition at r = is given by the Born distribution (2.4) which vanishes for p ab ^0. This 
implies that I (pab, t, v) vanishes for p a \ } — > and that the kernel is regular for p aq — > or 
Pqb —> 0. To make explicit the regularity of the kernel we use the splitting (for simplicity 
we replace p a b^P, Paq^Pi and p qb ^ p 2 ) 

p 1 , 1 , p - p\ - p2 



Pi P2 



1 1 
— + — + 

Pi P2 



Pi P2 



(2.7) 



where the last term is regular both for pi — > and p 2 ■ 



d T I(p, t) 



Using 



one finds 



r dn q 2 

J ~2k ~p~2 



d T I(p, t) 



dflq 
27T [p 2 [pi 



dpi 



I {pur) - I(p,t] 



0. The evolution equation becomes 
P-P1-P2 



2pi p2 



Hp, t ) 



(2.8) 



/' 

Jo 



\pi-p\ 



dpi 



I 



dQ q p - pi - P2 

27T pi p 2 



X dpi 
Pi 



^ , (2.9) 



J±I(pi, T )-I(p,T) 



+ 



J P Pi 



(2.10) 



\Pi ~ p\ VPi 

Using pi = r]p in the range < pi < p and pi = p/rj in the range p < pi < 1, one gets 
the final form of the evolution equation for the inclusive QQ distribution resumming all 
leading soft logarithms 

d T i( P , t)= C ^- (v-'iivp, t) - i(p, t)) + f-^- (iiv-'p, t) - i( P , t )) . (2.11) 
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The lower limit r)> p in the second integral ensures that the argument of I(p/rj, r) remains 
within the physical region p/rj < 1. 

In [3] the solution of (2.11) was considered only for p< 1. This corresponds to taking 
the e + e~ system in a boosted frame however it was not clear how to go back into a non- 
boosted frame. The reason for this limit is that it was argued that (2.11) reduces to the 
BFKL equation. Indeed, if one neglects the lower limit r] > p in the second integral of 
(2.11), the equation is formally the same as the BFKL equation 

d T T(p,r)= [ 1 -^L( v ^T(r ] p,T)-T(p,r)) + f -*L (T^ 1 p,r) - T(p,r)) , (2.12) 
Jo 1 - V Jo 1 - V 

where T(p, r) is the high energy elastic scattering amplitude in the impact parameter 
representation, p = b 2 is the square of the impact parameter and r = N c a s /irY (here a s 
is the fixed QCD coupling and Y the rapidity). 

Neglecting for p <C 1 the lower limit in the second integral of (2.11) one can use well 
known results in high energy physics such as the Pomeron intercept ap = %(0) N c a s /-ir 
where 

x (k) =2^(1) -^+ik)-^-ik), X (0)=41n2, (2.13) 

is the BFKL characteristic function. In this way one determines the asymptotic behaviour 
for large r of the inclusive distribution 

e x(o) r 

I(p,r) ^. (2.14) 

Neglecting the lower bound in the second integral in (2.11) implies that the inclusive 
distribution I(p, r) needs to be defined in the full range p > including the region p > 1 
outside the physical range. The extension to the non-physical region p > 1 seems not 
crucial since at high r the distribution rapidly decreases at large p so that one expects 
that such an extension should not give a significant contribution to the second integral 
of (2.11). A second argument supporting this expectation is obtained by observing that 
the QQ-equation is invariant if we change p^ap and by sending a — > one could neglect 
the lower bound in the second integral of (2.11). It was however pointed out by Gavin 
Salam that a cutoff in p actually gives even asymptotically an effect and generates [10] 
an additional factor r _1 in the large r behaviour in (2.14). This fact is discussed in the 
following. 

3. Comparison with BFKL equation 

Let us first recall the features of the solution to the BFKL equation which are relevant for 
the discussion of the solution of (2.11). The impact parameter b in the elastic scattering 
amplitude at high energy T(b 2 , r) is taken in the full range 6>0. In terms of the function 

o (x,r) =e x / 2 e- iln2T T(b 2 ,T), b 2 = e~ x , -cx)<x<cx), (3.1) 
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the BFKL equation (2.12) can be written in the form 

d T (f>o{x, t)= dy , — = -H ■ 4>q . (3.2) 

J-oo 2sinh 5 y — x 



2 sinh 2 |y — x| 

The kernel is translation invariant, hence diagonal in momentum space 

Jky _ e ikx 

2sinh±|y-x| ^ uvv " ^ 
E (k) = 4 log 2 -*(£;) = 14C(3) fe 2 + • 



/oo iny tux 
d y^~ TTl r = £o(£:)e fc 
-oo 2sinh 5 y — x 



The BFKL equation is then equivalent to a Euclidean Schroedinger equation (a diffusion 
equation) for a "free" particle with energy (dispersion relation) Eo(k) and the solution is 
given by a wave packet 

^e^C (k)e- TEo{k) , (3.4) 

where Co(k) is defined in terms of the initial condition. At large r the distribution is 
dominated by the small k region and consequently the solution is what we expect for a 
diffusion at large times, i.e. typically a Gaussian with standard deviation oc r 

CXD f fo -10 ) 2 1 /*°° 

00 (a, t) ^Co(O) P V-^- r J , C o (0) = / dx0 o (x,O). (3.5) 

with Z) = 28C(3). Translational invariance implies that the integral of the distribution is 
conserved 

d f°° 

— / dx(j)o(x,T) = 0, (3.6) 
dr J_ 00 

which is obviously the case also for the asymptotic expression (3.5). 

We come now to discuss the QQ-equation. As in (3.1), in the jet physics case we 
introduce the function 

0( X)T ) = e x/2 e~ 4ln2r I( P ,t), p = e~ x , < x < oo , (3.7) 

and (2.11) becomes 

d T 4>(x, r) = dy ^y^-^r) _ g ^ + ^ x/2) t)= _ h ^_ (3 _ g) 
Jo 2smh2|y-x| 

To discuss the solution of (3.8) we can still conveniently follow the analogy with the 
Schroedinger equation, but we have to take into account two different features: the exis- 
tence of a boundary at x = and a short range repulsive "potential" V(x) = 2 \og(l+e~ x / 2 ) 
which is responsible for the decrease in time of the integrated density 

J />00 poo 

— dx<f>{x,r) = -2 / dx\n(l + e~ x / 2 )<p{x,r) < 0. (3.9) 
dr Jo Jo 

(If we prefer to think in terms of diffusion processes, V is connected to a probability density 
for killing the particle). At large x (i.e. small p) the effect of the boundary and the potential 
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become (almost) negligible and the particle is essentially free with energy Eo(k) in (3.3). It 
will be natural to discuss the problem in terms of stationary waves, which should be given, 
at least for large x, by a shifted wave sin(/cx + S(k)), as is well-known from elementary 
wave mechanics. 

In order to explore this picture without technical complications, we shall first study, 
in the next section, a solvable model with and without boundary at x = 0. 

4. A simplified solvable model 

We consider the following evolution equation which bears some features similar to inclusive 
jet physics cases in (2.11) 

d T l(p, t) = f 1 dq r}- l l{r]p) + T dr)l(rr X p) . (4.1) 
Jo J p 

This is obtained from (2.11) by neglecting the singularity at rj = 1 and, correspondingly, 
neglecting the regularization due to the virtual correction. The crucial point of the dis- 
cussion here is the presence of the cutoff p < 1 which is reflected in the lower bound in the 
second integral. 

To discuss the effect of the cutoff p < 1 we consider the associated evolution equation 
(which bears some resemblance with the BFKL equation) for a function T(p, r) with p in 
the full range p > 

d T T(p, t) = f 1 dr, vT x T{j}p) + T dr? Tfa" V) • (4.2) 
Jo Jo 

Introducing as in (3.1) the function 

tpo(x,T) = e x/2 e~ iT T(p,t) , p = e~ x , -oo<x<oo, (4.3) 

we obtain (4 In 2 is replaced by 4 so that the integrated distribution of is conserved) 

f°° l 

d T Mx, t)= dy e-2\y-*\ (i; (y, T ) - fo(x, r)) = -Ho ■ Vo • (4.4) 



Due to translation invariance the evolution equation is diagonal in momentum space and 
the dispersion relation is now 

H ■ e ikx = £ (k) e ikx , £ (k) = - (4.5) 



(1 + 4 A: 2 ) 



thus giving the asymptotic solution (3.5) with D = 32. 
Consider now Eq. (4.1). We introduce the function 



${x, t) = e x/2 e- 4 T I(p, r) , p = e~ x , < x < oo , (4.6) 

and obtain 

r°° i 

d^(x,r) = / dy e~2\y~ x \ (^(y,r) - ip(x,r)) - 2e~^ ^{x,t) = —H • ip . (4.7) 
Jo 



-7- 



Since at large x the effect of boundary and potential V(x) = 2e~ x l 2 become negligible, the 
particle is essentially free with energy £o(k). At large x the solution should be given by a 
shifted wave sin(/cx + 6(h)). This model is actually solvable and this expectation turns out 
to be correct. 

To solve the model, observe that the kernel of (4.7) is easily identified with the Green's 
function of the operator — d/dx 2 with boundary condition at x = 

c> x V(0,t) = ±V(0,t). (4.8) 

This implies that the a complete set of eigenfunctions is given by 

. sin(A;x) + 2k cos(kx) . ., i , , A . 

u(x,k) = — i — ' 1 — - = sm(kx + S(k)), S(k) = tan -1 2k . (4.9) 
VI + 4k 2 

with the same eigenvalue £o(k) of (4.5). The solution is then 

*I>{x,t) = yj^ J dku( Xl k)C{k)e- T£ ^ k \ C{k) = yjl J dxu(x,k)tp(x,Q) . 

(4.10) 

For instance, by taking I(p, 0) = p (see (2.4)) onehasC(fe) = 8k/(l+4k 2 ) 3/2 . The behavior 
of the solution at large r is dominated by the small k behavior of the integrand. Because of 
the presence of the boundary the eigenfunctions vanish linearly at k = 0, also the amplitude 
C(k) vanishes linearly 

C(k) = c k (I + O (k 2 )) , (4.11) 

and then for large r 

with D = 32. The prefactor x+2 ensures that the asymptotic solution satisfies the bound- 
ary condition (4.8). The solution at large r has an additional factor {Dt)~ 1 w.r.t. the free 
evolution case in (3.5) (see also [10]). 



5. Solution of the jet equation 

The discussion on the solution of the complete problem for the inclusive off-shell gluon 
distribution I(p, r) will follow the same path as in the model case. We have to take into 
account the effect of potential V(x) and the presence of the border at x = 0. Since the 
potential is short range, for a function with support far from the boundary, the operator 
H in (3.8) is indistinguishable from the BFKL one in (3.2); it follows, like in ordinary 
quantum mechanics, that the continuous spectrum is unaffected. Also, since the potential 
is repulsive, we do not even have bound states. On the basis of what we have learned from 
the solvable model, we expect that the eigenfunctions u(x, k) are given by 

H ■ u(x, k) = E (k) u(x, k) 

(5.1) 

u(x, k) ~ sin(kx + 5(k)) , x>l, 
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with Eo(k) given by the BFKL dispersion relation (3.3). 

To find the expression for the phase shift S(k) and the eigenfunctions at finite x may 
seem problematic since, unlike the simplified model, there does not exist any simple bound- 
ary condition to be imposed to the eigenstates. Both the phase shift S(k) and the eigen- 
functions u(x, k) however can be determined by an expansion method which can be pushed 
to any order. Here we outline the method which will be reported in detail in Appendix A. 

One starts by defining the eigenfunction at large x as the shifted waves in (5.1). A 
first estimate of the phase shift is then obtained by requiring that the correction at large 
x is exponentially small. This also requires to determine the subasymptotic correction to 
the eigenfunction. Applying the operator H in (3.8) to a plane wave we get 

poo ik(y~x) _ i 

H e ikx = e ikx / — r dy - 2 log(l + e~ x/2 ) e ikx 

Jo 2s'mh±\x-y\ 

/oo ik(y—x) 1 poo —iky i 
— r -dy-e ikx / ^dy-2\og{l + e' x l 2 )e ikx 
.oo 2sinh±|x - y\ J x 2sinh±y 

= E (k) e ikx - r -A^ e -/ 2 + o {e~ x ) 

The leading term O [eT x l 2 ^ can be canceled by taking a superposition of e lkx and e~ lkx . 
The right superposition is precisely sin(/cx + So(k)) with So(k) = tan -1 2k the phase shift 
in the simple model (see (4.9)). To this accuracy the eigenfunction is given by 

uo(x, k) = sin(/cx + So(k)) , (5-2) 

with a correction of order e~ x . Of course, this is just the first step. To cancel the e~ x 
mismatch with the exact eigenfunction one should modify the phase shift and include a 
e -x correc tion in the wave- function with a coefficient which is determined by requiring that 
the residual correction is of order e~ Zx l 2 . This gives also the modified phase shift. The 
residual correction of order e~ 3x / 2 can be eliminated by repeating the procedure. 

In this way one sets up a recurrence procedure to determine both the phase shift and 
the eigenfunction, as it is shown in Appendix A. We introduce an expansion of the form 

u N (x,k) = Im |exp{-i(/cx + S N (k))} ^1 + ^c n e" n:E ^ j (5.3) 

where the coefficients {c n } and the phase shifts <5jv(A;) are determined in such a way that 
the error obeys the relation 



(H - E (k))u N (x,k) = O (e-^W*^ . 



(5.4) 



In Appendix A we show that the formula at finite can be extrapolated at N = oo and 
we get a closed formula for the phase shift 



X n^ ,, o, , 1 i \T(l-2ik)T 2 (l + ik) 

= ^ * k ( 41 ° g2 " 2C(3) k2 + ° ^ 



(5.5) 
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We also get a general expression for the approximate eigenfunctions at all orders which can 
be extrapolated to N — ► oo yielding the orthonormal eigenfunctions in terms of hypergeo- 
metric functions: 

u(x, k) = Im je^+W) 2 p x (\ -ik,\- ik; 1 - 2ik; e~ x ) } . (5.6) 

This representation is suited for approximating the large x region, but not for x^O corre- 
sponding to the physical 1 (the back-to-back dipole system for e + e in the center 
of mass). Notice indeed that the hyper geometric function diverges logarithmically for x^O 
but the phase shift is precisely tuned to cancel the divergence. To see this, we use one of 
Rummer's transformations for hyper geometric functions [11] to get a representation of the 
wave function which is manifestly regular at the origin. We can write Eq. 5.6 in the form 

u{x, k) = M(k) e~ ikx 2 F 1 (± + ik, \ + ik; 1; 1 - e~ x ) . (5.7) 



with Af(k) = y/kn tanh/c7r. Notice that, in spite of the factor e tkx , the eigenfunctions are 
real 3 . This representation is well suited for small x corresponding to p close to one 

u(x, k) = N{k) {l + {\-k 2 )x + (x 2 )) . (5.8) 

In conclusion, for a given initial condition l(p,0), the QQ-multiplicity is given by 

<f>(x, r) = y - J dk u(x, k) C{k) e~ T Bo(fe) , 

C(k) = ^J^ dxu{x,k)e% l(e~ x ,0). 



(5.9) 



The dispersion relation Eo(k) is given by the expression (3.3) as in BFKL case. The initial 
condition at r = is given in (2.4). For instance, taking I(p, 0) = ^fpp^ with /i>0we have 

™-W^£f" (5-10) 

The eigenfunction expansion can be reduced to the well-known Mehler-Fock transform 
[11, 12] by expressing u(x,k) in terms of Legendre functions. In terms of the original 
density I(p, r) one then finds 



I( P ,T)=fdki(k)P_ h+tk (?f) 



(5-11) 

I(k) = ktanhirk J dz I \j^,0) P_i +ik (z) 

In Fig. 1 we show, as function of r, the values of <j>[x, r) at x = (i.e. p = 1) corresponding 
to e + e~ in the center of mass system. Here we show also the corresponding QQ-distribution 
I(p=1,t) (we have considered a simplified initial condition). 

In Fig. 2 we show <fi(x, r) as function of x for values of r which are experimentally 
accessible. The value at x = decreases by increasing r as seen in Fig. 1. 

3 This alternate presentation of the eigenfunctions was pointed out to us by V. A. Fateev. 
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From (5.9) one has that the solution for large r is dominated by the region around k = 0. 
The phase shift is a simple analytic function of k which vanishes at k = 0. Therefore also the 
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coefficient vanishes linearly with k. Since the product u(x,k) C(k) vanishes quadratically 
for k — > 0, at large r the distribution is given by the derivative of a Gaussian, which 
implies the additional power r _1 as expected (see also [10]). Fig. 3 shows that the Gaussian 
behaviour starts to set up already at moderate values of r. 

In Appendix B we describe a numerical approach to obtain the solution of the QQ 
evolution including the large r region, which, in practice, is faster than using the exact 
solution. At large r the shape is very accurately fitted by 

(x + x {t)) 2 ' 



(x,t) 



x + x (t) 
A - exp 



(r 



to) 2 



2D(t - tq) 



(5.12) 



where D = 28 £(3) and x$ is positive and weakly r-dependent; this feature is beyond the 
leading saddle point approximation. As in the model previously discussed (see (4.12)) this 
expression vanishes at a negative value (outside the physical region). 



x 10 




Figure 3: Long time evolution of 4>(x,t) obtained from (5.9) with initial condition l(p,0) = \p 



6. Conclusions 

The QQ and BFKL equations (see (2.11) and (2.12)) are formally similar but originate 
from multi-soft gluon distributions in two completely different phase space regions. This 
fact produces important differences. First of all while in the BFKL equation (see (3.2)) the 
variable x ranges on the entire axis, in the QQ case (see (3.8)) it is bounded to be positive. 

Moreover, due to the fact that qu (the soft gluon transverse momentum) is the ar- 
gument of the QCD running coupling for each emitted gluon, the evolution parameter r 
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is given by (2.2) in the QQ case and by a s Y in the BFKL case. Since in BFKL case 
a s remains fixed, r is large at high energy and then it is physically relevant to study the 
asymptotic solution of the BFKL equation. In the QQ case one has that r involves the 
running coupling at the hard scale and then increasing Q one has that r increases slowly. 
So the physically accessible values of r are not large. 

A further difference between the two cases is in the physically relevant values for the 
variable x. In the QQ case for e + e~ annihilation in the center of mass one has x = 0. In 
the high energy case one is interested in small impact parameter (short distance) which 
corresponds to large x. 

We have given the solution for the QQ equation which can be used in the physical 
range (finite r and small x). The inclusive distribution I(p,r) giving the multiplicity for a 
QQ pair is 

!(P, r) = ^j o °° dk u(x, k) C(k) \ p= , (6.1) 

with C(k) given in (5.9), r given in (2.2), the BFKL characteristic function in (2.13), 
the eigenfunction u(x,k) given in Eq.(5.6-5.7) and l(p,0) the Born distribution given in 
(2.4). Here 6 is the angle of the dipole of partons emitting the QQ system. In e + e~ anni- 
hilation in the center of mass the dipole is given by the back-to-back quark-antiquark and 
one has 9 = it. In performing these integrations one may find difficulties due to the oscil- 
lating form of the eigenfunctions. An alternative is to use directly the evolution equation 
(3.8) which can actually be easily studied numerically. We describe in the Appendix B the 
method we used and the check of the analytical results. 
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A. Analytic solution of the spectral problem 



We start by constructing a systematic expansion as outlined in Sec. 5. To this aim, it is 
technically more convenient to use the original representation Eq. (2.11). Let us define 



— A+n+ik 

— P* 



We have 



H 4> n = / — : dr] (j) n + / — dr] (j> r , 

JO V J p I — T] 

= H n {k) <p n + V n (p, k) <p n 



where 



.(*) = f E ( 



Pn(k)= I > ' ( rT +n+ik ^ - V m ) d V 



E 



m + n + i + ik m + l 



m>0 

ip(l) -il)(\+ik + n) 



and 



: l rt~2- n - ik - 1 
v n (p,k) = / dr] 

1 — 7/ 



1 1 \ ^ / p \+m-n-ik pm +l 



= y I L_ )-y _ 



, +m-n + ik m+l J z -' \\+m — n + ik m + l 

m>0 \ ^ / m>0 \ ^ 

The first sum amounts to ^(1) — ^ Q — n — ik) which can be added to \i n to give 

Xn{k) = p n + v n = 2^(1) -tp(^ + ik + n) — ip (^ — ik — n) 

Note that xo is just Lipatov's function Eq. (3.3), up to an additive constant. The basic 
relation we are to use is then the following: 

„m+n+3/2+ik „m+l 

HMp,k) = xn(k)Mp,k) + £ - — n E t— r (A- 1 ) 

^— ' m+l ^ ± + m — n — ik 

rn>0 rn>0 ^ 

We now define 

$ N (x, fc) = Im J e"* 5 °- MP, k) \ , (co = 1) (A.2) 

I n=0 J 

which fails to be an eigenfunction for H because of the various terms originating from 
Eq. (A.l). The coefficients c n and the phase shift 5 will be determined by requiring that 
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the difference 



R N (p,k) = {H- X o{k))$N 

n>l 



, L » m+n+3/2+ifc 

^ ^ m+ 1 ^ 

n>Ora>0 

-vy -; — ^ — P m+i 

be of the highest order possible O (p( N+3 ^ 2 ). 

For the simplest possibility, N = 0, we only have S to be fixed in order to cancel the 
leading term in Rq which is — ~ This gives the same phase shift of the simplified 
model of Sec. 4, namely 5o(k) = tan -1 2k. 

Consider now N = 1. We have 

R 1 (p, k) = (1 + Cl ( X i - xo)) P 3/Wfc + O (p 5/2 ) 
^ | - ifc — ifc^) 

Now we fix ci in order to cancel the (p 3 / 2+lfc ) term and the phase shift to cancel the 
O (p) term: 

l + 2ik . n , , 2/c 

Ci = , 01 = tan 2/e + tan — . 

In general, for AT even, we shall determine ci, C2, . . . , cjv/2 by canceling the terms of order 
(p 3 / 2 ) , C (/? 5 / 2 ) , . . . ,0 ( /5 ( iV + 1 )/ 2 ) and the coefficients c N/2+1 ,c N/2+2 , ■ ■ ■ , c N from the 
terms of order O (p 2 ) , . . . , O (p 1+N ^ 2 ) . The remaining term, the leading one of order 
O (p), is canceled by fixing 8. Similarly, we can proceed for N odd. The coefficients c n are 
given by the recurrence relation 

£™=l C n-m/™ co = 1) (A3) 



where we have used the relation x n (0) — Xn(&) = ^Y^=l (2j + 1 + 2ifc) 1 . From the first 
few c n it is natural to conjecture 4 that 

{\ + ik) n 2 

Cn ~ n\(l + 2ik) n (A - 4) 

where as usual (a) n = T(a + n)/T(a), the Pochhammer symbol. The expression for the 
higher coefficients requires the solution of a linear system and is more cumbersome. The 



The proof of the conjecture has been found independently by G. E. Andrews and V. A. Fateev, private 
communication. 
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phase shift at the lowest orders is given by 



Si(k) 


= tan~ 


1 2k + tan 


63(h) 


= tan~ 


1 2k - tan 


S 5 (k) 


= tan~ 


1 2k - tan 


S 7 (k) 


= tan~ 


1 2k - tan 



2k 

T 

_ 1 2/c 1 2/c 

k + 2 tan 1 h tan 1 — 

3 7 

1 2/c _ 1 2/c 1 2/c 

fc + tan" 1 h 2 tan -1 h tan 1 — 

3 5 11 

1 /c 1 2/c -1 2/c 

« — tan — h tan h 2 tan — 

z O o 

1 2/c 1 2/c 

+ 2 tan" 1 — -f-tan" 1 — 
7 15 

A; 2k 

5g(k) = tan _1 2k — tan -1 k — tan -1 — + tan" 1 — 

1 2/c „ 1 2/c -1 2/c 1 2/c 

+ tan" 1 — + 2 tan" 1 — + 2 tan" 1 — + tan" 1 — . 
5 7 9 19 

The regularity of the pattern allows to extrapolate to any N, namely 

n k 2n 2k 

^2n+i(fc) = tan" 1 2k — > tan" 1 h > tan" 1 

^— ' m ^— ' 2m+l 

m=l m=± 

Ei 2/c 1 2/c 
tan" 1 + tan" 1 

9™ 4-1 8n 4- 



(A.5) 



(A.6) 



2m+l 8n + 3' 

m=ra+l 



a similar relation holding for 821% ■ By taking the limit n — ► 00, we get the following result 
for the phase shift 

5{ k ) = 4 log 2 k + f; (-r 2(22r " 2 1 r ) + C i (2r+1) fc2r+1 ^ A - 7 ) 

r=l 

which is precisely the expansion of Eq. (5.5). In the limit N — > 00 we get the eigenfunction 
representation given in (5.6). 

At this stage a conservative statement is that this equation represents just an asymp- 
totic approximation to the solution, owing to the fact that the higher order coefficients 
( c tv/2+i> •••) c n) grow quite rapidly with N. If we examine the finite sums $tv of Eq. (A. 2) 
near x = we find a typical behaviour of asymptotic approximations, namely they oscillate 
around the function u(x, k) of Eq. (5.6). 

Our conjecture is that the solution in terms of hyper geometric functions is exact, but 
an interchange of integration and power series expansion is not allowed. This obstruction 
may be due to the fact that the integral operator H is unbounded. From the agreement 
with the numerical calculations of Appendix B and the overall consistency of the whole 
picture, we argue that this problem can probably be ignored, but a rigorous proof would 
be welcome 5 . 



5 An exact solution for a class of integral equations, including the one considered here, and related to 
Tuck's equation [13] has been found in the meanwhile by V.A.Fateev (private communication, see [9]). 
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We can cure the divergence of the finite approximation by applying Pade approxi- 
mants, which effectively regulate the oscillating terms. Notice that the hyper geometric 
function diverges logarithmically for x — ► but the phase shift is precisely tuned to cancel 
the divergence. To see this, we use one of Rummer's transformations for hyper geometric 
functions [11] to get the representation of the wave function given in (5.7), which is mani- 
festly regular at the origin. From this representation we can extract the intercept at p = 1 
as a function of k to be compared to the numerical value obtained with the methods of 
Appendix B. This agreement gives a further support to the consistency of the analytic 
solution. 

B. Numerical study 

The evolution of the density I(p, r) can be easily studied numerically in the representation 
given by Eq. (3.8). By introducing a normalization interval < x < L and a discrete 
grid xo = 0, x\,X2, ... ,xn = L with uniform spacing a = L/N, the integral operator Hq 
can be approximated by a finite N x N matrix to various degrees of accuracy. A simple 
choice is to adopt a trapezoidal rule to represent the integral; the singular nature of the 
kernel does not raise serious difficulties. One can choose a tolerance parameter e, typically 
1CT 8 — 10~ 12 and put all matrix elements less than e to zero, in such a way that the 
matrix can conveniently be implemented as a sparse matrix. We checked that an initial 
wave- function with support near x = L/2 evolves according to the BFKL equation until it 
becomes sensitive to the boundary. 

The evolution in r was studied using the routines provided by matlab. By varying L 
and a we determined that L = 300 and a = 0.1 are acceptable, in the sense that finer grids 
and larger L do not vary the results appreciably. Taking <j>(x, 0) = e~^ x one can compute 
4>(x, r), for < r < 100 in a few seconds on a desktop computer. The snapshots in Figs. 1-3 
are obtained this way, but they are indistiguishable with what one can obtain using the 
solution (5.9). 

At large r the shape is very accurately reproduced by 

a( \ ~ a x + x o(r) f (x + x (t)) 2 \ 

* A v^r exp V m^)\ ■ ( } 

Using D = 28^(3) and fitting the other parameters, we find that a coincides with | to 
within 0.2%. Notice that this is not a universal property of the solutions. If the initial 
shape has support at large x then its evolution is effectively described by BFKL at least 
until boundary effects set in. 

The program can also investigate the eigenvalues and eigenfunctions. Clearly, since 
we work in a finite interval, the spectrum is discrete with spacing w n/L. From this 
preliminary numerical investigations it emerged the nature of the eigenfunctions as phase- 
shifted standing waves, and we were able to check that i) the spectrum is indeed given by 
Eq. (2.13), and ii) the phase shift vanishes for small k. 
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